echo on
% script file per confrontare vari algoritmi per il calcolo dei valori
% del polinomio (x-1)^6
%
p = [ 1 -6 15 -20 15 -6 1];
x = linspace(0.995,1.005,1000);
pause
% Algoritmo di horner
y = horn(p,x);
pause
z = (x-1).^6;
pause
% Base delle potenze
w = pol(p,x);
pause
plot(x,y,'.r',x,z,'-g',x,w,'.b')
